Quest for rare events in mesoscopic disordered metals 
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The study reports on the first large statistics numerical experiment searching for rare eigenstates 
of anomalously high amplitudes in three-dimensional diffusive metallic conductors. Only a small 
fraction of a huge number of investigated eigenfunctions generates the far asymptotic tail of their 
amplitude distribution function. The relevance of the relationship between disorder and spectral 
averaging, as well as of the quantum transport properties of the investigated mesoscopic samples, 
for the numerical exploration of eigenstate statistics is divulged. The quest provides exact results 
to serve as a reference point in understanding the limits of approximations employed in different 
analytical predictions, and thereby the physics (quantum vs semiclassical) behind large deviations 
from the universal predictions of random matrix theory. 
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Despite years of extensive studies, initiated by the 
seminal work of Anderson, 1 the problem of quantum 
particle in a random potential ("disorder") still poses 
new challenges. The biggest shift in approach and in- 
tuition came with the advent of mesoscopic quantum 
physics 2 — from the critical phenomena like description 3 
of the localization-delocalization (LD) transition, to re- 
alization that complete understanding of quantum inter- 
ference effects, generating LD (zero-temperature) tran- 
sition for strong enough disorder, requires to study the 
full distribution functions 4 of relevant physical quantities 
(such as conductance, local density of states, current re- 
laxation times, etc. 5 ) in finite-size disordered electronic 
systems. Even in the diffusive metallic conductors, the 
tails of such distribution functions show large deviations 
from the ubiquitous Gaussian distributions, expected in 
the limit of infinite dimensionless zero-temperature con- 
ductance g = G/(e 2 /Trh). These asymptotic tails are 
putative precursors of the developing Anderson localiza- 
tion, occurring at 3 g ~ 1. 

A conjecture about unusual eigenstates, being micro- 
scopically responsible for the asymptotic tails 4 of various 
distribution functions related to transport, was put for- 
ward early in the development of the mesoscopic pro- 
gram. 5 However, it is only recently 6 that eigenstate 
statistics have come into the focus of mesoscopic com- 
munity. Similar investigation in the guise of quantum 
chaos 2,7 started a decade earlier (leading to the concept 
of scarring 8 in quantum chaotic wave functions). Thus, 
the notion of "prelocalized" states has emerged from the 
studies of quantum disordered systems 5 ' 9,10 of finite g. 
In three-dimensional (3D) conductors these states have 
sharp amplitude peaks on the top of a homogeneous back- 
ground. Their appearance, even in good metals, has been 
viewed as a sign of an incipient localization. 5 However, 
the non-universality of the prelocalized states tempers 
this view, 11 nonetheless, pointing out to the importance 
of non-scmiclassical effects 11 in systems whose typical 
transport properties are properly described by semiclas- 
sical theories [including the perturbative quantum cor- 
rections 12 ~ 0(g^ 1 )]. Here I provide an insight into 



the weirdness of such states: Fig. 1 plots the amplitude 
spikes of an anomalously rare prelocalized state, to be 
contrasted with the amplitudes of an ordinary extended 
state shown in the same Figure. The study of prelocalized 
states is not only revealing mathematical peculiarities of 
the eigenproblem of random Hamiltonians, but is relevant 
for various quantum transport experiments, such as e.g., 
tunneling experiments on quantum dots where coupling 
to external leads depends sensitively on the local prop- 
erties of wave functions. 13 By exploiting the correspon- 
dence between the Schrodinger and Maxwell equations 
in microwave cavities, it has become possible to probe 
directly the microscopic structure of quantum chaotic or 
disordered wave functions. 14 

This study presents a numerical result for the statis- 
tics of eigenfunction "intensities" ^^(r)! 2 in closed 3D 
mesoscopic conductors which are diffusive L 3> I (L and 
£ being the size of the system and elastic mean free path, 
respectively), metallic g 3> 1, and weakly disordered 
kpl ^> 1 (hp is the Fermi wave vector). The statisti- 
cal properties of eigenstates are described by a disorder- 
averaged distribution function 10 

fE {t) = T^)N ^ - \**(r)\ a V)5(E - E a )\ , (1) 

on N discrete points r inside a sample of volume V. Here 
p(E) = {^ a 5{E — E a )) is the mean level density at 
energy E, and (. . .) denotes disorder-averaging. Finite- 
size disordered samples are modeled by a tight-binding 
Hamiltonian 

# = ^e m |m)(m| + ^ *«nn|m)(n|, (2) 

m <m,n) 

with nearest-neighbor hopping i mn = 1 (unit of energy) 
between s-orbitals (r|m) = ip(r — m) located on sites m 
of a simple cubic lattice of size L = 12a (a being the lat- 
tice spacing). Periodic boundary conditions are chosen 
in all directions. The disorder is simulated by taking the 
potential energy e m to be a uniformly distributed ran- 
dom variable over the interval [— W/2, W/2]. This is a 
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FIG. 1. An example of eigenstates in the band center 
of a delocalized phase in the Anderson model of W = 5. 
The average conductance at half filling is g m 12, entailing 
anomalous rarity of the "prelocalized" states. For plotting 
of the eigenfunction values in 3D, the sites m of 12 3 lattice 
are mapped onto {1,...,1728} in a lexicographic order, i.e., 
m = (m x ,m y , m z ) ^ 144(m x — 1) + 12(m H — 1) + m z . 



standard Anderson model of localization, 1 where weak 
disorder is introduced by W — 5 (the model can also 
be viewed as a discretized version of the problem of a 
single particle in a continuous random potential). The 
eigenproblem of (2) is solved exactly by numerical diag- 
onalization. A small energy window AE — 0.007 is po- 
sitioned around E = 0, picking up 3 states in each 

of the 30000 conductors with different impurity configu- 
rations. Altogether, 43432 states have been collected for 
the evaluation of /e=o (*) as a histogram of intensities at 
all points inside the sample (N = 12 3 ). This large en- 
semble makes it possible to obtain a well-defined far tail 
(Fig. 2), stemming from the amplitudes of prelocalized 
states (the far tail has been investigated by numerical 
simulation only in low-dimensional systems 15 ' 16 where, 
strictly speaking, all states are localized even for weak 
disorder, 3 albeit with an exponential localization length 
in 2D). 

In both generic quantum chaotic and quantum dis- 
ordered systems eigenstates are characterized solely by 
their energy (the only constant of motion) instead of a 
set of quantum numbers. Their classical counterparts 
are nonintegrable — in the entire phase space for "hard" 
chaos. 7 Quantum chaos implies semiclassical approxima- 
tion 7 (ft — > 0), and is less efficient 14,16 in localizing 
chaotic wave functions through scars than localization 
by quantum disorder, 1 due to interference effects in the 
diffusive motion. Since eigenstates and eigenvalues can- 
not be obtained analytically, one usually resorts to some 
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FIG. 2. Statistics of eigenfunctions in the band center of 
the Anderson model with diagonal disorder W = 5 on a sim- 
ple cubic lattice 12 3 . The numerically calculated distribution 
f E=0 (t) is fitted with: leading order of the Fyodorov-Mirlin 
correction, f PT (t)[l + 0.04(3/2 - 3t + t 2 /2)], to PT distri- 
bution (valid in the perturbative region); and exp-log-cube, 
0.0256 exp(— 0.239 hi 3 £), at large deviations from /pt(£) in 
the non-perturbative region — similar fit is achieved with 
exp-log-cube, 0.045 exp[-0. 19 ln 3 (i/0.72)], akin to / N lsm(*). 



statistical treatment. The unifying concepts in this pur- 
suit come from approaches like random matrix theory 7 
(RMT) and its justification (note that matrix elements 
of Hamiltonians of disordered solids are spatially depen- 
dent and therefore do not satisfy standard assumptions 
of RMT) through supersymmetric nonlinear cr-model 12 
(NLSM). RMT gives universal predictions for the level 
and eigenstate statistics, which are applicable in the limit 
g — > oo. Its answer for fsit), in time-reversal invariant 
systems, is given by the Porter-Thomas 7 (PT) distribu- 
tion function 



frr(t) 



1 



: exp(-t/2). 



(3) 



Within RMT, x f , Q (r) is a Gaussian random variable, 
thereby leading to x 2 distribution of intensities. The 
NLSM allows one to study the deviations from the univer- 
sality limit in the finite-size weakly disordered (fep^ 3> 1) 
samples, which turn out to be determined by the prop- 
erties of classical diffusion. 6 The distribution function 
fE=o{t) is contrasted with fpr(t) in Fig. 2. For t > 5, a 
deviation in the tail starts to develop, eventually becom- 
ing a few orders of magnitude greater probability than 
RMT prediction for the appearance of states with high 
amplitude splashes. However, the disorder is so weak 
that large-i limit tail is shorter, and composed of smaller 
values of /^(i), than some quantum chaotic tails at- 
tributed to scarring. 8 ' 16 Small deviations from PT distri- 
bution are taken into account through perturbative (weak 
localization) correction 18 
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FIG. 3. Comparison of disorder average (over the impu- 
rity ensemble) and combined disorder and spectral averaging 
[over small interval A_E used to broadened S(E — E a ) from 
Eq. (1) into a box function] in computing /e(£): (a) histogram 
of the number of eigenstates which exhibit the largest ampli- 
tude splashes determining the very end of the far tail of the 
distribution of |\& a (r)| 2 — both histograms are defined by only 
1.4% of 43342 or 1.86 • 10 6 eigenstates collected in the inter- 
val AE = 0.007 or AE = 0.3 around E = 0, respectively; 
(b) perturbative region and FM fit (note a small discrepancy 
between the data and the fit, which is increasing upon ap- 
proaching the upper boundary W ~ 6 of the semiclassical 
transport regime); (c) far tail region and exp- log-cube fit. 



derived by Fyodorov and Mirlin (FM) 18 for t <C k -1 / 2 . 
Here k (~ L/gi in 3D) has the meaning of a classical 
(time-integrated) return probability for a diffusive par- 
ticle. 6 In the region of considerable deviations, where 
decays much slower than /pt(*)> one has to use 
the non-pcrturbative predictions for the asymptotics of 
/#(£) in mesoscopic metals. 12 They employ either NLSM 
techniques 6 



/nlsm(*) ~ exp 



1 
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t > k^ 1 



(5) 



or "direct optimal fluctuation method" of Rcf. 1 1 (where 
only the leading term of exp-log-cube asymptotics is eval- 
uated explicitly 11 ) 



/dof(*) ~ exp [-Cdof k F £ In 3 1] 



(6) 



aiming to describe the short-scale (and non-semiclassical) 
structure of the solutions of Schrodingcr equation in a 
white-noise disorder. The intermediate region of ampli- 
tudes (Fig. 2) is covered in the NLSM formalism by 



nt 2 



i<t<i. 



(7) 

where a correction in the exponent is large compared to 
unity, but small compared to the leading RMT term. 10 

Numerical techniques provide the means to solve ex- 
actly the (discretized) Schrodinger equation for a quan- 
tum particle in an arbitrary strong random potential, i.e., 
for all transport regimes — from semiclassical to strongly 
localized. 19 Such solutions are a good reference point 
from which one can gauge analytical predictions and their 
"neglect" of those effects which are hard to deal within 
a specific formalism. Analogous numerical investigations 
°f m quantum chaotic systems have typically dealt 
with a single set of eigenstates 16 since impurity ensemble 
is absent in standard clean examples of quantum chaos 
(like e.g., quantum billiards 14 ). Even in quantum dis- 
ordered systems, the computation of has involved 
only relatively small ensembles of disordered conductors, 
or combined disorder and spectral averages (the disorder- 
averaging is supposed to improve the statistics, while 
claiming that individual eigenstate would show more or 
less the same behavior 16 ). Thus, the natural questions 
arise: Is an enormous ensemble of disordered conductors 
indispensable to obtain complete in a good metal? 
Can spectral averaging [i.e., wide broadening of S(E—E a ) 
in Eq. (1) into a box function of the width AE] be used to 
facilitate this quest? Both questions demand an empiri- 
cal answer, which is readily obtained from Fig. 3. Only a 
small fraction of all investigated eigenfunctions exhibits 
the highest observed amplitudes [panel (a) in Fig. 3]. 
Therefore, one has to search for special configurations 
of disorder where quantum interference effects can lead 
to large wave function inhomogencitics. They arc anoma- 
lously rare in weakly disordered conductors, i.e., the tail 
of /b(^) is not generated by a large number of more or 
less equal- weight impurity configurations. 11 The effect of 
spectral averaging is not a priori obvious. 20 Namely, the 
short scale physics to be discussed below can correlate 
eigenstates with neighboring eigenvalues and therefore 
produce artifacts of the numerical procedure employed 
to get fE(t) in a finite-size system with discrete spec- 
trum, when the original definition (1) in terms of delta 
functions is abandoned. Nonetheless, the comparison of 
fE(t) from Fig. 2 and the one obtained from the same im- 
purity ensemble but with AE = 0.3 (where 62 states are 
collected, on average, from each sample) suggests that 
combined impurity ensemble and reasonable spectral av- 
eraging might be the most efficient way to obtain a long 
and smooth far tail (Fig. 3). 

Mesoscopic transport properties of finite-size samples 
can be fully delineated from the "measurement" on a 
computer. This offers a simple way of comparing the 
above formulas to the observation in real systems — 
instead of trying to deduce the functional form of /#(£) 
phcnomcnologically (which is usually inconclusive be- 
cause of the possibility to fit reasonably-well different 
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functions), one can use the "measured" values, adjust the 
free parameters, and, if the analytical expressions work, 
explore the viability of physics behind the raw numbers. 
The quantitative understanding of transport is also pre- 
requisite for another reason. The lattice size sets a lower 
limit on the disorder strength which ensures the diffusive 
transport regime (L»l). When disorder is too strong, 
semiclassical parameters become nonsensical (£ < a), al- 
though the conductor might still be far away (j » 1) 
from the LD transition point. 3 ' 21 On the other hand, only 
for strong enough disorder the far tail is reasonably long 
(Fig. 2) to allow the extraction of reliable fitting param- 
eters. The interplay of these three limits leaves a narrow 
range 4 < W < 6 of possible disorder strengths in this 
study. 

The exact zero-temperature disorder-averaged conduc- 
tance is calculated from the two-probe Landauer formula 
g(E F ) = Tr[t{E F )tf(E F )] w 12 (at E F = 0). 21 The 
Fermi wave vector k F w 2.8/a, averaged over the Fermi 
surface E F = of a simple cubic lattice, serves as a coun- 
terpart of k F used in theoretical simplifications which 
assume a Fermi sphere. For W = 5 disorder g is still 
dominated by the semiclassical effects. 21 Thus, I use the 
Bloch-Boltzmann formalism (applicable when k F £ 3> 1), 
in Born approximation for the scattering on a single im- 
purity, to obtain 21 I w 1.4a (k F £ w 4). 

Small deviations (both in the body and in the tail) 
of the numerical fE=o(t) from the PT distribution arc 
well-explained by /fm(£) [panel (b) in Fig. 3]. However, 
k ~ 0.08 extracted from the fit is much larger than the 
diffusive and universal 18 (i.e., independent of the details 
of disorder), K diff = (2/gn 2 ) £ q exp(-X>q 2 r)/q 2 L 2 = 
0.018. The exact value of the sum over the diffusion 
modes (q being quantized by periodic boundary condi- 
tions) is evaluated numerically where exponential cut- 
off provides the necessary ultraviolet regularization at 
|q| ~ £ _1 (V = v f t/3 is a classical diffusion constant). 
The divergence of the sum in 3D necessitates to go be- 
yond the diffusion approximation of the NLSM, imply- 
ing that short-range spatial correlations become impor- 
tant. This is to be contrasted with the long-range cor- 
relations (described by the soft modes of NLSM) in low- 
dimensional systems where Kdiff is convergent parameter- 
free number. Therefore, to reproduce the fitted n, the 
other relevant ("ballistic" and non- universal) contribu- 
tion Kbaii ~ (k F £) , arising from quantum dynamics on 
time scales shorter than elastic mean free time 6,22 r, has 
to be included n = Kdis + «baii ~ (k F £)~ 2 + (kpi)^ 1 . 
The intermediate region is poorly described by the for- 
mula (7). 

The beginning of the far tail t > k _1 ~ 12 is lo- 
cated approximately using the fitted k from the pertur- 
bative region. The exp- log-cube formulas (5), (6) im- 
ply that this region should be explained by an exponen- 
tial function of a cubic polynomial, C p exp(— C 3 In 3 t — 
C2 In 2 1 — C\ lnt). The fit, however, cannot give the ac- 
curate values for all four parameters simultaneously. On 



the other hand, the fit including just the cubic term 
C p exp(— Csln 3 ^) works well (Fig. 2, 3), giving C 3 = 
0.239 and C p = 0.0256. The NLSM universal result goes 
beyond predicting only the prcfactor of the leading log- 
cube term by including the lower powers of Int. This 
form also fits the data reasonably well with k = 1.39 and 
log-cube prefactor Cnlsm ~ 0.19 [which is close to the 
prefactor 1/4/t = 0.18 in /nlsm(^)]- Nevertheless, there 
are two major puzzles with /nlsm(^) fitting the numeri- 
cal data: (1) inasmuch as k is a characteristic quantity of 
classical diffusive dynamics, determining both perturba- 
tive and non-perturbative corrections (within the frame- 
work of NLSM) to RMT picture, 6 the different values, re- 
quired here for different intervals of amplitudes in Fig. 2, 
are hard to reconcile (regardless of the dynamical origin 
of k); (2) since Kdiff ~ (k F £)~ 2 , the log-cube prefactor 
has the form 1/4k <~ (k F £) 2 . However, the investigation 
of different ensembles for a range of I (oc W~ 2 ) hints to- 
ward the prefactor (C 3 or Cnlsm) being close to linear 
function of k F i as is the case with 23 /dof(^) prediction 
(an unexpected sublinear prefactor behavior is seen also 
in 2D Anderson model 15 ). 

Thus, the exp-log-cube formulas could account for the 
far tail of statistics of eigenstates in the Anderson model. 
However, their expected parameters are inadequate, sug- 
gesting that short scale effects below £, which are to 
be treated fully quantum mechanically, 11 ' 22 are essen- 
tial (even in the perturbative region and for large sam- 
ples). They lead to strong dependence of eigenstate 
statistics in 3D on microscopic details of a random po- 
tential. 6 It remains to be seen if the log-cube prefactor 
can be explained by examining the short-scale physics 
(thereby going beyond standard diffusive NLSM correc- 
tions, which are semiclassical in nature and universal) — 
Anderson model can hardly satisfy some key assumptions 
of the "direct optimal fluctuation method", 11 while bal- 
listic NLSM, 9 as yet, has only reproduced the diffusive 
NLSM result. 11 
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